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Abstract 

We predict contributions of higher i^-resonances to the radiative rare decays b 
S7, in the framework of the QCD sum rules on the hght cone, which combines the 
traditional QCD sum rule method with the description of final state mesons in terms 
of the light-cone wave functions of increasing twist. Our calculations are restricted 
to the leading twist-two operators for i^*(892) and to the asymptotic wave function 
for the other if **-mesons. 

Using experimental data on the semileptonic r — > K**Vr decays, we extract the 
corresponding decay constants for vector and axial-vector iC**-mesons. Based on our 
estimate of the transition form factor F^*^^'^'^\q) = 0.32 ±0.06 and the if|(1430) de- 
cay constant /i^|(i43o) = (160 ±20) MeV , we find good agreement with experimental 
results. The two largest fractions of the inclusive b ^ branching ratio are found 
to be (10.0 ± 4.0)% for B K*{892)'j and (5.0 ± 2.0)% for B Er*(1430)7 decays. 
We also compare our results with the existing theoretical predictions. 



*e-niail : safir@mail.desy.de 
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1 Introduction 



The flavour-changing neutral current (FCNC) -B-decays, involving the 6-quark transi- 
tion 6 — i> (s, (i) + 7 and b (s, d) + l^l^ {I = e, fi, r, z/) provide a crucial testing grounds 
for the standard model at the quantum level, since such transitions are forbidden in the 
Born approximation. Hence, these rare -B-decays are characterized by their high sensitivity 
to New physics. 

In the standard model, the short distance contribution to rare B-decays is dominated 
by the top quark, and long distance contributions by form factors. Precise measurements of 
these transition will not only provide a good estimate of the top quark mass and the CKM 
matrix elements Vtd, Vts, Vtb, but also of the hadronic properties of 5-mesons, namely form 
factors which in turn would provide a good knowledge of the corresponding dynamics and 
more hint for the non-perturbative regime of QCD. 

The experimental searches for the FCNC B-decays have already provided first evidents. 
The recent experimental observations of the rare decay mode B — > K*^ have been deter- 
mined by CLEO and more recently also by BABAR [§ and BELLE 



and 



(4.55 ± 0.70 ± 0.34) X 10-^ [1] 
5r(fi° ^ ir*°7) = <; (4.39 ± 0.41 ± 0.27) X 10-^ § (1) 

(4.96 ± 0.67 ± 0.45) X 10-5 § 



and also of the inclusive rate ^ ^ 

Br{B X,7) = (3.22 ± 0.40) x 10"^ (3) 

However, the first observation of the rare S-decay to the orbitally excited strange mesons 
has been reported only by CLEO [Q]. With an integrated luminosity of 9.2 fh~^ of e^e~ 
data corresponding to 9.7 x 10^ BB mesons pairs, they measure the B i^'g (1430)7 decay 
with a branching fraction of 

Br{B i^*(1430)7) = (1.66+°:| ± 0.13) x 10"^ (4) 

These important experimental measurements provides a crucial challenge to the theory. 
Many theoretical approaches have been employed to predict the exclusive B fC* (892)7 
decay rate (for a review see |^ and references therein). On the other hand less attention 
has been devoted to rare radiative 5-decays to excited strange mesons [p|-p!l|]. Most of 
these theoretical approaches rely on non-relativistic quark models HQET |10| and 



relativistic model |]TT|. However there is a large spread between different results, due to a 
different treatment of the long distance effects. 

The aim of this paper is to estimate the contribution of higher ii"-resonances to the 
radiative rare decays b 57. While the effective Hamiltonian approach allows to include 
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short distances QCD corrections at scales /i^ ~ m^, estimates of the hadronic matrix 
elements of the relevant operators necessarily require some non-perturbative technique. 
We use here QCD sum rule |jl2|, which has been proved to be a powerful tool for such 
purposes. 

In this paper, we use an elaborate version of the QCD sum rule, known as QCD sum rule 
on the light cone which has been developed originally for light quark systems in [|T3],Q. 
This method was demonstrated to be very suitable for exclusive process. Moreover, it 
is technically much simpler than the standard QCD sum rule approach, where the basic 
correlator between the vacuum and a light meson state is expressed by universal quantities 
called light-cone wave functions for that particular light meson. These wave functions 
represent distributions of the light-cone momentum of the constituents, and can be classified 
by their twist defined as the difference between the canonical dimension and the Lorentz 
spin of the corresponding operator. Their asymptotic form is fixed by perturbative QCD, 
while the non-asymptotic effects at lower momentum scales can be estimated from QCD 



sum rules for two point correlators of appropriate currents |15, 16 



The principal theoretical result of this paper is a sum rule for the electromagnetic 
penguin form factor appearing in the decays B — > K**'~f, where K** can be pseudoscalar, 
vector, scalar, axial-vector, tensor or pseudo-tensor. We derive the corresponding sum rule. 

The paper is organized as follows. In section 2, we review the effective Hamiltonian for 
radiative -B-decays, dominated by the electromagnetic penguins. In section 3, we define 
the form factors, which govern the exclusive rare -B-decays to orbitally excited i^'**-mesons. 
In section 4, we describe the method of the QCD sum rules on the light cone and derive 
our general sum rule for the form factors of the electromagnetic penguin operators entering 
in the decays B V{A)'y and B T{Ta)^ where V{A) and T{Ta) stand respectively 
for vector (axial-vector) and tensor (pseudo-tensor) i^**-meson. The wave functions and 
the decay constants which are needed to evaluate the sum rule are discussed in section 5, 
and the numerical estimates of the form factors and the predicted Branching fractions of 
B K**'y are presented in section 6. Finally, section 7 contains a summary and some 
concluding remarks. 



2 Effective Hamiltonian 

At the quark level, the rare semileptonic decay b s'y can be described in terms of 
the effective Hamiltonian obtained by integrating out the top quark and bosons: 

H^fj = -A^V*Vt, J2 C^{^^m^^) . (5) 

* i=l 

where Vij are the corresponding CKM matrix elements and is the Fermi coupling 
constant. The operator basis is defined as follows |]T^,[T^: 

O2 = {sLalfihlfi) {cl^Yclo) 
q=u,d,s,c,b 
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q=u,d,s,c,b 

O5 = (sLalfibla) ^ {qRfil^qRp) 
q=u,d,s,c,b 

Oe = (sLali^bLfs) ^ {qRpl^qRa) 
q=u,d,s,c,b 

g 

'^7 = T? 9 ^" '^Mi' ("^bR + TTLsL) ha F^" 

iO TT^ 



where L/R = (1 =1= 75)72, a^^ = f [7^,7jy] and = {ps —p)fj. is the four-momenta transfer 
of the photon. T", a = 1....8 are the generators of QCD, and a,P are SU{3) color indices. 
Here F^'^ and G°'t^'^ denote the electromagnetic and chromomagnetic field strength tensor, 
respectively. The Wilson coefficients Cj(/i) are evaluated perturbatively at the W scale 
and then evaluated down to the renormalization scale fi ~ rrib by the renormalisation 



group equation. Actually, they are known at the next-to-leading order . 



3 Matrix elements 

Our purpose in this section is to write the specific matrix elements for B K**'~f 
transitions, where K** denotes the higher ii'-resonances summarized in table |I|. At the 
tree level only the magnetic moment operator O7 contributes to the transition amplitudes. 
Thus, we are interested only in the following matrix elements 



(/r*|e,.gn5), {K^QUlB) (7) 

where = sa^^b and 9^^ = sa^ul^b are respectively the tensor and pseudo-tensor 
currents. For a generic radiative decay B — >• F7, where F = P,V, S, A, T and Ta stands for 
pseudoscalar, vector, scalar, axial- vector, tensor and pseudo-tensor respectively, one defines 
a transition form factor F[^{q^) as: 



{P{Sm\Q%{Q,M^\B{pB)) = (8) 

(P(S)(p)|0,.(eJjgn5M) = ^ ^l^^J^^^ {q\pB+p), - (ml - ml^,^)q,] (9) 

iy{A){p,e)\Q,^{QlM'^\B{pB)) = te,,,^ e*'' p'^ 2F^^^\q') (10) 
{ViA)ip,e)\Qlie,^)q-^\BipB)) = F^^^\q') {e;(m| - ml^^^) - {e^ps) {pb+p),] 

+ F^^''\q^Wps) \q, - , {PB+p)An) 

I ^B ^V{A) I 



where ^^^^^^(O) = F^'^^\o), and: 



4 



(r(T^)(p,6)|e,.(eyg1i?(pB)) = ^e,^,^ e*''- p'^ 2F^^^^\q') (12) 

\ niB ^ mB \ 



X e 



P% Pb 



a0 2 



{{ml - m|(T^^)) Qf, - q'^{pB+p)^} (13) 



where rriF, €^ and e^'^ are respectively the generic mass of the i^'**-meson, the polarization 
of the vector (axial-vector) and the tensor (pseudo-tensor) final meson. 

With the above definition, the exclusive decay widths are given by : 

T(B^P{Sh) = (14) 

T(B ^ V(Ay,) = Onl\CAm,)\^FrW ( 1 " f 1 + (15) 

\ niB I \ niB 

r(B^T(T,H, = i„=|C,(m,)PFr™(0)'^fl-^-''' 

6 ™T{Ta) \ ^B 

l + ^^l (16) 

with ( = "^2tt^ l^b^fsP- A good quantity to test the model dependence of the form factors 
for the exclusive decay is the ratio of the exclusive-to-inclusive radiative decay branching 
ratio |n]]: 

^^^^ BR{B-^Xsi) ' ^ ^ (1 -m2/m52)='(l + m2/mfe2) ^ ' 

^ ^ BR{B ^ T{Ta) 7) _ ^^^^(0)2 ' ^t(t^)) + ^t(Ta), 

^^^^^ ~ BR{B^Xsl) ~ ^ml^^TA) (1-^2/771^2)3(1 + ^2/771^2) 

where 777ir = niF/mB- With this normalization, one eliminates the uncertainties from the 
CKM factor Vt^VJ* and the short distance Wilson coefficient Cj{mb). Thus, we are left 
in eqs.(|l^ and (|18D with unknown form factors F/'(0), which we will derive in the next 
section. 



4 QCD sum rules on the light cone 

For definiteness, we show here and in the next section the derivation of the sum rule 
for the transition form factor -^/^(O) as defined in eqs.(p!0|) and (|12|), governing the radiative 
decays B K**'y. 
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Meson 






Mass (MeV) 


K 


0~ 




497.67 ±0.03 


K*{892) 


1- 


I'Si 


896.1 ±0.26 


/s:i(1270) 


1+ 


l^Pi/l^Pi 


1273 ± 7 




1+ 


l^Pi/l^Pi 


1402 ± 7 


K*luiO) 


1- 


23^1 


1414 ± 15 


K*{U30) 


0+ 




1412 ±6 


7^2* (1430) 


2+ 




1432 ± 1.3 


K{im) 


0' 




^ 1460 


7^2(1580) 


2- 




^ 1580 


ii'i(1630) 


?• 




1629 ± 7 


/Tifieso) 


1+ 


[2^P,] 


1650 ± 50 


ir*(1680) 


1" 




1717 ±27 


i^2(1770) 


2- 




1773 ± 8 


X3*(1780) 


3- 




1776 ± 7 


i^2(1820) 


2- 




1816 ± 13 


i^'(1830) 


Q- 




^ 1830 


fs:o*(1950) 


0+ 


[2'Po] 


1950 ± 10 ± 20 


ii'2*(1980) 


2+ 


2'P2 


1973 ± 8 ± 25 


ii'*(2045) 


4+ 




2045 ± 9 


7^2(2250) 


2- 




2247 ± 17 


ir3(2320) 


3+ 




2324 ± 24 


X5*(2380) 


5- 




2382 ± 14 ± 19 


7^4(2500) 


4- 




2490 ± 20 


ir(3100) 


?• 




^ 3100 



Table 1: Spectrum of excited neutral K-meson taken from /|^/, except for the quantities in 
brackets. The quantum numbers n, L, s and J represent respectively the radial excitation, 
the orbital angular momentum, the sum of the two quark spins in the meson and the total 
meson spin. The parity is defined as P = (—1)'^+^. 
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The starting point of our sum rule is to consider the correlation function 



I / dx e*^^ < K**{p,e)\T {i,{x)a^,{a,,^,)q''h{x)m 275^^(0)} |0 > (19) 



which we will study in details for the vector (axial-vector) and tensor (pseudo-tensor) K**- 
mesons. Hereafter we use as a generic notation for the field of the light quark. Our 
approach is very close to the calculation of the B — > K{K*) form factor in [p] ,^. 

4.1 B ^ ^(^) form factors 

To derive the sum rule, we start with the following matrix element of the time-ordered 
product of two currents between the vacuum state and the vector (axial-vector) i^**-meson 
at momentum p: 



I J dx e'"^ < V{A){p,e)\T {ij{x)a,,ia,,j5)q''bix)m zj^i^iO)} \0 > = i ^e^^'^'^g^'^ 

X T((p + g)2) (20) 

at = and at Euclidean ml — {{p + q)^) of order several GeV^. The hadronic representa- 
tion of (PO) is obtained by inserting a complete set of states including the i?-meson ground 



state, higher resonances and the non-resonant states with i?-meson quantum numbers 



KA). p ^ ^2^ _ {V{A){p, eMx)a,A^,,^5)q''b{x)\B) (i?|&(0) ^75^(0)|Q) 

te^,,^e qp l[[p + q) ) - ^|_(p + g)2 

^ ^ {V{A){p,e)\s{x)a,,{a,,^r:,)q''b{x)\h) (/i|&(0) ^75V^(0)|0) 
^ m| - (p + qy 

(21) 

Writing down the dispersion relation in {p + q^, we can separate the contribution of the 
i?-meson as the pole contribution to the invariant function T{{p + q)"^), i.e. 



nip + qf) = ^ 2 f J,2 + 22 

where the dots stand for contributions of higher-mass resonances and the continuum. The 
i?-meson decay constant is defined in the usual way, 

< 0|^7M75&|i?(p) >=wJb, (23) 
rriB and rrih are the 5-meson and the 6-quark masses, respectively. 
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The possibility to calculate the correlator ([20| ) in the region of large space-like momenta 
(p + g)^ < is based on the expansion of the T-product of the currents in (|20| ) near the 
light-cone = 0. The leading contribution to the operator product expansion arises from 
the contraction of the 6-quark operator in ( pO] ) to the free 6-quark propagator < 0|66|0 >. 
The light quark operators are left uncontracted. Diagrammatically, this contribution is 
depicted in Fig.(l.a). The formal expression is easily obtained from (|20|): 



dx e 



iqx 



q 



(2vr) 



< V{A){p,e)\T {ij{x)a^,{a^,^5){mb+ h) ^5^(0)} |0 > (24) 



where we made use of the following representation of the free propagator S^ix): 



< 0\T {b{x)b{0)} |0 > 



X] 



-ipx 



mb+ j) 
p2 — ml 



da 



la 



(25) 



In general, the expression ( P^ is expressed through matrix elements of non-local op- 
erators, sandwiched in between the vacuum and the meson state. These matrix elements 
define the light-cone meson wave functions. The first term in (|2^) is given by ||2^ : 



<0|V>(0)a^,(a^,75)V'(^)im)(P,e) > = ^(e(,") p. - p^) /^(^) 



du e-™P^ 



(26) 



Likewise, 



<0|^(0)7^(7^75)^(x)|\/(A)(p,e) > = fv(A)mv(A) du e-'^^^ x 

Jo 

f (X) (e^^^a;) 



e(%) 
{p x) 



{p x) 



M I 2\ 



v(A)mv(A) 



X / d«e-^"^^^?it(^)(«,/i') 

'0 



\,ViA){u, ^M^) 

(27) 



f28) 



The functions 4>±y{A){u, fi'^) and 0||y(A)(M, /x^) give the leading-twist distributions in 
the fraction of total momentum carried by the quark in transversely and longitudinally 
polarized mesons respectively. The functions (/^^^^^^(m, /i^) and g±^Y(^A)('^^ 1^'^) twist-3 
wave functions. The normalization is chosen in such a way that for all all four distributions 

/ = (p±,v{A), 9^ly{A) and ^^y(^), we have: 



du f{u) = 1 



(29) 
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In the matrix elements of non-local operators on the l.h.s of (^)-(|28|) the separations 
are assumed to be light-like, i.e. = 0. Regularization of ultraviolet divergences that arise 
in the process of the extraction of the leading — behavior produces a non-trivial scale 
dependence of the wave functions, which can be found by renormalization group methods. 
The scale in (^) is fixed by the actual light cone separation, /i^ ~ ~ ml — {p — qY- 

Putting eqs.(P^)-(PBD together, we obtain: 



du 

""6 



1 



ml + uuml^^^^ 



u{p + g)2 



X 



(30) 



+ 



du- 



mr 



2 2 

^V(A) 



■^:^;^^:^^^V{A)fviA)9^y^A) 



4-/o " '{ml + uuml^^^ 



u 



where u = 1 — u. The expression in (^) has the form of a dispersion integral in {p + q)"^, 
which can be made explicit by introducing the squared mass of the intermediate state 
s = ml/u + umy^^^ as the integration variable, instead of the Feynman parameter u.The 
basic assumption of the QCD sum rule approach is that the contribution of the i?-meson cor- 
responds in this dispersion integral to the contribution of intermediate states with squared 
masses smaller than a certain threshold sq (duality interval). Making the Borel transfor- 
mation l/(s — (p + q)"^) — * exp{—s/t) and equating the result to the 5-meson contribution 
in (|22|), we arrive at the sum rule 



+ 



du —exp 
u 



u^m; 
t u 



+ m 



V{A)j 



e 



- 2 

So umy,y. 

u ^ ' 



'mbfv{A)4>±y{A){u, fl =t)+ Umv(A)fv{A)9ly^A)^U, /i = t) 
+ ut _ 

(31) 



m 



2 2 



mv(A)fv{A)9±y(A)i^^f^'^ = 



Aut 



which should be satisfied for values of the Borel parameter t of order several GeV'^. We 
see that Borel transformation removes arbitrary polynomials in {p + qY and suppresses the 
contributions from excited and continuum states exponentially relative to the ground-state 
contribution. 



4.2 B T{Ta) form factors 

The aim of this subsection is to calculate the transition form factor f'^'^'^^^ as defined 
in ([I^), governing the radiative decays B —>■ T{TA)'y. We consider first the correlation 
function: 



F,{q,p) = 1 dx e''^%T{TA){p,e)\T{4,{x)a^,^,{<J,,)q'b{x)m Z75^(0)} |0) (32) 
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The hadronic representation of (^) is obtained in the same way as for (20). Finally, 
one gets 



F,{g,p) = I /rfxe— ^-^f. (35) 



F,{q,p) = e:/^F{q^ {p + qf) p, + e*J^H{q\ {p + qf) q, 
^ niB rnB 

+ el^—G{q\{p + qf) (33) 
^ rriB 

where q denotes the four-momentum transfer. From now on, we shall concentrate on the 
invariant amplitude F which is physically more interesting than the amplitude H and G. 
For F one can write a general dispersion relation in the momentum squared (p + q)"^ of the 
5-meson: 

(34) 

— [P + q) 

From eqs.(p5D and (0), one can obtain the leading contribution of Fig.(|^), and equals 

da 

X ^T(T^)(p,e)|T|v^(x)cr^,75K-)g"(m, + ^^ ^75^(0)| |0 

From now on, we will restrict ourself just to the helicity A = for the tensor and 
axial-tensor mesons, due to the fact that this polarization brings mainly the leading twist-2 
distribution. The other polarizations contribute in (^) with more higher twist distribu- 
tions, which will not be considered here. 

Let us consider the first term in (|35D . The matrix element of the nonlocal operator is 
given by |T5[: 

<T(rA)(p,e)A=o|V^(x)a^.75K.)g"75^(0)|0>=-2p^/T(T^) ! rfw e*"''^0T(T^)(M, /i^) (36) 

^0 

where (f)T{TA){'^^ 1^"^) is the tensor (pseudo-tensor) ii'**-meson light-cone wave function of 
the leading twist-2, normalized as in (pQ]). Substituting the matrix elements (^) into (p5D, 
and integrating over x and the auxiliary parameter a we find the following expression for 
the coefficient p^ in (|3^ ) that is for the invariant amplitude F: 



F{q'^, {p + qf) = -fT(TA)T^b I du 

Jo 



{p + q uY — ml 



(37) 



After using the Borel transformation of (pTQ and equating to (|3J), the transition form 
factor f'[^'^^^ is obtained 



^T(T^)(0) 



B JA 



^ du 
— exp 



u 



m 



B 



ut 



0t(Ta)(m,/x^) 



(3^ 
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where A = (m^/so) and friTA) denotes the tensor (pseudo-tensor) decay constant. 

In addition to the quark-antiquark wave functions considered above there are in principle 
also contributions from multi-particle wave functions. The most important corrections 
of this type are expected to arise from quark-gluon operators in the operator product 
expansion of (|2^) and (|32D . A typical diagram where the gluon is emitted from the heavy 
quark is shown in Fig. (lb). These contributions bring twist-3 contributions, which we do 
not take into account in this paper. 

Principal inputs in these sum rules are the vector (axial-vector) and the tensor (pseudo- 
tensor) meson wave functions and the meson decay constants, which encode non-trivial 
information about the dynamics at large distances, and which we are going to present now. 



5 Wave functions of the ir**-mesons 

Here, we collect the formulae for the light-cone wave functions of the i^**-mesons and 
specify the parameters. It is important to note that the asymptotic form of these functions 
and the scale dependence are given by perturbative QCD The twist-2 wave function 



,K*{892)iu, fj, ) is expressed as an expansion in Gegenbauer polynomials 



K*{892){U, fi^) = Qu{l - U) 



1 + ai(/i2)e + a^il?) ( - J) + a,{i^) " ^ ) + 



(39) 



where we have introduced the shorthand notation ^ = 2u — 1. The non-perturbative 
effects are contained in the coefficients a„. In the leading logarithmic accuracy, they are 
multiplicatively renormalizable and have the following scale dependence: 



aM = <^n{^^o)[^^y' (40) 



."s(Aio), 

where h = {ll/3)Nc — (2/3)r2j is the LO coefficient of the QCD beta function, Nc and Uf 
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being the number of colors and active flavours, respectively. The anomalous dimensions 
turn out to be 



ji+i 



In 



(41) 



where Cp = {N^ — l)/2A^c- The non-asymptotic coefficient in ([39| ) can be estimated from 
two-point sum rules for the moments J vP' (p ±,k* (892) {u, ^'^)du at low n, or taken from 
experiment. The non-perturbative information encoded in the quark and gluon condensates 
is thereby transmuted into the long-distance properties of the wave function. For the 
numerical results, we use the following estimate at /i^ = 1 GeV^ [p!5|,p2 



af = 0.75, af(^^^) = -2, af'^'"'^ = 0.75 (42) 

As mentioned above, the correct normalization point in the sum rules is the order of the 
typical Borel parameter, which for i?-meson decays is of order /i^ ~ — ml ~ 6 GeV^. 
Using Ag^ = 0.225 GeV with (0) and (H), one gets for the evoluted coefficients 



a 



_ft:*(892) 



0.57, 



a 



A'* (892) 



-1.35, 



a 



K*(892) 



0.46 



(43) 



The wave functions corresponding to the above parameterization are shown in Fig.(^). 
As a general effect of the rescaling, the wave functions become somewhat wider and closer 
to the asymptotic expression (p±{u) = 6u{l — u). In the case of the iir*-meson, the wave 
function becomes both wider and more symmetric, so the scaling violation affects the region 
M < 1/2 only. Happily, the integration over u in the sum rule in (|31|) is restricted to the 
interval of rather large momentum fractions, carried by the light quark involved in the 
electromagnetic penguin operator. For realistic values of parameters, one has for the K*- 
meson u > 0.61 — 0.7. According to the analysis in the i^'*-wave function turns out 
to be in this region very close to the asymptotic expression. 



The conformal expansion of the wave functions g^f^j^ 



(892) 



somewhat more involved and can be obtained using the approach |14 



and fl'^_^^.(892) 



U, /i J IS 

It should be noted 

that these wave functions contain contributions coming from both operators of twist-2 and 
twist-3. The twist-2 contributions to the transverse wave functions g^l^x*{892) x*(892) 



receive a Wandzura-Wilczek type contribution I^S] that can be expressed in terms of the 
leading-twist longitudinal distribution function 0||,_ft'*(892): 



9±,K*(892)y'^) 



^^5±,if*(892)l 



U 



dw 



4>\\,K* (892) (w) 



W 



1 — w 



dw— 
1 

dw 



|,A-*(892)lWj 
W 

0||,i^*(892)(w 
W 



(44) 



In our evaluation we will not consider the twist-3 contributions to 9±^k*{892)('^^ 1^'^) 
9±^K*{892)('^^ which can be written in terms of three-particle antiquark-gluon-quark 



12 



1.5 



0.5 



-0.5^ 





/ ^ — 
/ / 

/ / 


X 




- 










- 












/ r 










/ r 










/ r 










/ / ' 










/ / 1 












0.2 


0.4 


0.6 


0.8 


1 



u 



Figure 2: The leading twist-wave function (t)±_{u) from the sum rule analysis in [TSj, for 
the K*-meson. Solid and Dashed lines correspond respectively to the K*-wave function at 
the scale /i^ = 5.8 and 1 GeV^. The Thick one represents the asymptotic wave function 
= 6m(1 — u). 



wave functions of transversely polarized vector meson |T5|,p^. Following we write 
these functions 



9 1_^K* (892) i'^) 



:i + e 



(45) 



Concerning the other i^'**-meson wave functions the situation is more complicated. It 
should be noted that the higher resonances are notoriously difficult to handle using QCD 
sum rules, because it is not possible to get rid in the sum rule of the contributions of the 
lower-lying states with the same quantum numbers. One possibility to solve this problem, 
as it is shown in Fig.(0) is to use the asymptotic wave function for i^**-mesons. Since 
the light quark involved in the electromagnetic penguin operator for B K** transitions, 
carries a large momentum fractions u > 0.62 — 0.72 in the sum rule (P3), then according to 



ref. the corresponding K**-wave function in this region is very close to the asymptotic 
one. 

Thus, restricting our calculations for the ii'**-wave functions to the asymptotic mode, 
we do not consider the higher twist contributions 



u 



6 u (1 — u) 



(46) 
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1 777 GpV 




0.2 GeV 


a(5) 
^^QCD 






137 


Gp 


1.16637 10"^ 


tb 


1.61 ps 


tr 


0.291 ps 



Table 2: Default values of the input parameters used in our numerical calculations. 

<t>T(TA){u) = 6m(1-m) 

^V(A)(«) = (47) 
^V(A)(«) = 

Clearly, the asymptotic wave functions should be normalized to the characteristic scale of 
the process under consideration. Usually, this brings non-asymptotic effects into play which 
change the shape of the above wave functions, but preserve their normalization to unity. 
Unfortunately being of non-perturbative origin, these effects are difficult to evaluate. 



6 The decay constants 



The extraction of the 5-meson decay constant fs from QCD sum rule has received a 
lot of attention in the past [p5| -|33[|, using for the two-point correlator of 675M-currents. 

The sum rule in yields ~ 120 — 160 MeV , which is lower than the value preferred 
at present, typically ~ 180 MeV (see |^,|32| for a review). In fact this smaller value 



is an artifact of neglecting the radiative corrections, which are numerically large. Recently 
calculated perturbative three-loop QCD corrections are incorporated into the sum rule, 
leading to fs = 197 ± 23 MeV in ref. @ and ~ 206 ± 20 MeV in ref. |5|. Hence, we 
observe that the effect of the 0{as) correction to the value of fs is sizeable. 

We note that what is actually calculated in ( ^I]) and (|38D are sum rules for the convoluted 
product of two amplitudes, fsFflO). From refs. and [^, it is important to notice that 
F^*{0) and increase either with higher operator contributions or with 0{as) corrections. 
However, one can hope that going beyond the asymptotic distribution for Ff{0) and fs, 
the ratio of the corresponding sum rules will be stable. Since we are discarding the 0{as) 
radiative corrections in our sum rule and restricting our calculation to the leading twist-2 
for B i^* (892)7 decay and to the asymptotic distribution for B K**'^ decays, we will 
take a reasonable value of fs ~ 160 MeV for *^^^^^(0) and ~ 120 MeV for "(0), 
as a guess of our estimate on the 5-decay constant in our approach. 



For the K* (892)-meson decay constants, we adopt the same choice as in ref. . Their 
numerical values determined from QCD sum rules are |1T2|,|15|: 



fK' =210 MeV 



(4J 



In the case of i^'**-mesons, the situation is more complicated since their decay constants. 
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is:*(1270) Ki{UOO) K*{U10) K*{U30) fs:i(1650) ir*(1680) 



UMeV) 



1+ 
122 



1+ 
91 



1- 
86 



0+ 
79 



1+ 
86 



1- 
86 



Table 3: Central values of the pseudoscalar, vector, scalar and axial vector K** -meson 
decay constants 



are not known, and we are not aware of any estimates of their numerical values. One 
should extract them experimentally or rely on some models. Using the semileptonic r — > 
(P, y, S, A)ut- experimental data |2^, it is possible to extract the corresponding meson decay 
constants. These quantities contribute to the transition rates for pseudoscalar, vector, 
scalar and axial vector emission, and we find from straightforward calculations |]36[ 



/p(5) 

fv(A) 



'Br 



exp 



in rrir 



m 



V(A) 



mt 



\ 



Br 



exp 



in m 



V{A) 



m 



V{A) 



mt 



(49) 



(50) 



where t^- is the r-lifetime, the quantities fp{s) and fv{A) are respectively the pseudoscalar 
(scalar) and the vector (axial-vector) meson decay constants. Using the central values of 
^''^T-^{pv s A)vr lis^^*^ 1131' obtain in table ^ the corresponding decay constants. 

Concerning the tensor and the pseudo-tensor meson decay constants, it is difficult to 
extract them from the r decay data, due to the impossibility to define them from the 
vacuum-to-tensor (pseudo-tensor) matrix element. A somewhat different approach can be 
used to obtain these missing decay constants, either using the SU(3) symmetry breaking 
or from other processes. 

The /^gll^SO) decay constant can be estimated from the experimental data on {B — >■ 
(1430)7) reported recently by CLEO [0]. Assuming a stability of our sum rule around 
t = 20 GeV"^ and t = 25 GeV'^ , we plot in Figs.(||) and (H) the constraint on /x2*(i43o)- 
The prediction of the i^2(1580) decay constant is somehow unknown and neither SU(3) 
symmetry breaking nor other processes are helpful. For that we assume a stability of our 
sum rule around t = 20 GeV"^ and t = 25 GeV"^ , and we plot in Figs.(H) and (H) the 
corresponding fractions of the inclusive b ^ s'j branching ratio (|1^) as function of the 
decay constant. 



7 Numerical results 



Having the necessary formulae at hand we next explain our choice of the relevant 
parameters. Let us first consider the parameters which characterize the 5-meson channel 
that is niB = 5.27 GeV, fs, rrih and sq. The value of the pole 6-quark mass rrih and the 
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continuum threshold Sq, are varied in the hmits: 

nib = 4.6 - 4.8 GeV, 

So = 33 -35 Ge1/^ (51) 

Before giving numerical predictions on the form factors Ff{0) we must first determine 
the range of values for the Borel parameter t for which the sum rules ( |HT| ) and (|3^) can 
be expected to yield reliable results. Usually, the lower limit of this range is determined 
by the requirements on the terms proportional to n > 1 to remain subdominant. In 
(pT|) this concerns in particular the twist-3 wave functions, which increase rapidly at small 
t. The upper limit of the allowed interval in t is determined by demanding the higher 
resonance and continuum contribution not to grow too large. We have checked numerically 
for 5 < t < 10 GeV^ that the resulting value of ^^^^""(O) is practically independent of 
the Borel parameter, as it is shown in Figs.(0) and (^). 

For other form factors involving B V{A)'~f transitions, we have used the asymptotic 
wave function (^61) , which induces good stability for F^^'^\o) form factors, in the same 



Borel range. 

The stability plots for the form factors as functions of the Borel parameter are shown in 
Figs.(|7|-p!8D. The stability is in all cases good, and the variations of the continuum threshold 
within the limit (|5TD induce uncertainties in the values of form factors within less than 22% 
for F^^^\ Using the region in the Borel parameter 5 GeV^ < t < 10 GeV^, we extract 
the values of the vector and the axial- vector form factors 



^i^*(892)^Q^ = 0.32 ±0.06 

^Ki{i270)^Q^ = 0.14 ±0.03 

^Ki{i4oo)^Q^ = 0.098 ±0.02 

^i^'(i4io)^Q^ = 0.094 ±0.02 

^i^i(i650)^Q^ = 0.091 ±0.02 

^i^*{i680)^Q^ = 0.091 ±0.02 (52) 

Using the sum rule (|38|), and the experimental data on B ^ ^2 (1430)7 0) "we show in 
Figs.(§) and (|) the experimental constraint on the (1430)-decay constant. The shaded 
area represents the uncertainty on the 6-mass. Assuming a stability of our sum rule (^) 
at t = 20 GeV^ and 25 GeV^, we obtain 

//^|(i430)..o = (160±20)Mey (53) 

Using the central value of (p3D, it turns out that the tensor form factor F^^''^^^^\o), as it 
is shown in Figs.(|l9D and (|3)? leads to 

^i^j(i430),=o(Q) = 0.19 ± 0.04 (54) 

In Figs.(^) and (^, we have plotted the normalized branching ratio of the 7^^2(1580)- 
meson as function of its decay constant. The shaded area depicts mainly the mfe-related 
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uncertainty in our sum rule (^) at t = 20 GeV^ and t = 25 GeV^. It turns out from 
Figs.(^) and (|^) that the uncertainty on the corresponding normahzed branching ratio 
increases rapidly with the decay constant. 

The given errors in (|52D and (|5^ ) mainly come from the uncertainties in the 6-quark 
mass and the dependence on the Borel parameter. This leads for all form factors except 
, to an uncertainty of order 20%. We observe that for all decays the average value 
of the momentum fraction u under the integral in ( ^I]) and (^) is < m >~ 0.6 — 0.7. It 
should be noted that errors given here reflect uncertainties in the input parameters, but do 
not include a possible theoretical uncertainty of the method itself. 

Using our form factors and (p^, we estimate the corresponding ratio R = r{B 
K**^)/T{B Xg'f) for various _ft'**-mesons. We find a substantial fraction (6.0 — 14)% 
of the inclusive b ^ s'y branching ratio going into the K*{892) channel, and (3.0 — 7.0)% 



going into the Jr|(1430) channel. Both of these predictions are in good agreement with the 
recent experimental result of Rk*(892) ~ (17.59 — 11.63)% and Rk*{u30) ~ (4.23 — 2.98)%, 
presented in (P, (|D and d). 

In order to make comparison of our results with previous calculations, we have tabulated 
our results together with results of 0, 0, and [Tl| in table ^. As far as we know. 



these four works are the only ones that have dealt with radiative rare B decays into higher 
i^'-resonances. Concerning B f^* (892)7 and B —>■ (1430)7, our predictions agree 



with the values quoted in [jlOl and [|T^]. Note that our result are in general in much better 
agreement with |T0| than and 

With the exception of the i^* (892)7 and i^^^ (1430)7 channels, no other exclusive radia- 
tive processes have been identified so far. Thus, it is difficult to make more quantitative 
statements on the other i^'**-decay channels, especially when the decay in consideration is 
totally dominated by the long distance effects. 
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Figure 3: Experimental constraint on /xj(i43o) at t = 20GeV^ and Sq = 35GeV 
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Figure 4: Experimental constraint on fK2{U30) at t = 25GeV^ and Sq = 35GeV^. 
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Figure 5: Stability plots for /x2(i580) at t = 20GeV'^ and sq = 35GeV^ 



10- 



Rk2{1580) [%] 




oL 



0.08 0.1 0.12 0.14 0.16 0.18 0.2 

fK2(i58o){GeV) 

Figure 6: Stability plots for /a-2(i580) at t = 25GeV'^ and Sq = 35GeV^. 
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Figure 7: Stability plots for the sum rule in of the Borel parameters for B fT* (892)7 
at Sq = 35 GeV^. Solid, Dashed and Thick lines correspond respectively to mi,=4-6 GeV, 
4.7 GeV and 4.8 GeV. 
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Figure 8: Stability plots for the sum rule in (31) of the Borel parameters for B — > K*(892)7 
at Sq = 33 GeV^. Solid, Dashed and Thick lines correspond respectively to mi,=4.6 GeV, 
4.7 GeV and 4.8 GeV. 
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Figure 9: Stability plots for the sum rule in of the Borel parameters for B 1270)7 
at Sq = 35 GeV^. Solid, Dashed and Thick lines correspond respectively to mi,=4-6 GeV, 
4.7 GeV and 4.8 GeV. 
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Figure 10: Stability plots for the sum rule in (3J_) of the Borel parameters for B — > 
1270)7 "^0 = 33 GeV^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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Figure 11: Stability plots for the sum rule in ( ply of the Borel parameters for B 

1400)7 '^^ "^0 = 35 GeV'^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV. 
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Figure 12: Stability plots for the sum rule in ( pl\ ) of the Borel parameters for B — > 
1400)7 '^t Sq = 33 GeV^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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Figure 13: Stability plots for the sum rule in of the Borel parameters for B — > 

1410)7 "So = 35 GeV"^ . Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV. 
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Figure 14: Stability plots for the sum rule in of the Borel parameters for B — > 
1410)7 "5o = 33 GeV^ . Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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Figure 15: Stability plots for the sum rule in ( ply of the Borel parameters for B 

1650)7 So = 35 GeV'^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV. 
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Figure 16: Stability plots for the sum rule in ( pl\ ) of the Borel parameters for B — > 
1650)7 '^t Sq = 33 GeV^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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Figure 17: Stability plots for the sum rule in ( ply of the Borel parameters for B 

1680)7 '^^ ■So = 35 GeV'^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV. 
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Figure 18: Stability plots for the sum rule in of the Borel parameters for B — > 
if* (1680)7 "5o = 33 GeV^ . Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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Figure 19: Stability plots for the sum rule in ( p^ j of the Borel parameters for B 
i^a (1430)7 "So = 35 GeV"^ . Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV. 
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Figure 20: Stability plots for the sum rule in ( ^ ) of the Borel parameters for B — > 
iir2 (1430)7 "5o = 33 GeV^. Solid, Dashed and Thick lines correspond respectively to 
mb=4.6 GeV, 4.7 GeV and 4.8 GeV 
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ref. |T0[ ref. g ref. 



ref. iO 



K 







forb. 








(892) 




10.0 ± 4.0 


16.8 ±6.4 


3.5- 


12.2 


4.5 


15 ± 3 


TV'* / 1 /I on\ 

K (14dUj 


u 






forb. 






Ki[121Q) 


1 + 


2.0 ± 0.8 


4.3 ± 1.6 


4.5- 


10.1 


forb./6.0 


1.5 ± 0.5 


iil(i4UU j 




O.y ± 0.4 


2.1 ±0.9 


6.0- 


13.0 


iorb./6.0 


2.6 ± 0.6 






O.U ± z.U 


6.2 ±2.9 


17.3 - 


37.1 


O.U 


0. / ± i.z 


K*(1680) 




0.7 ±0.3 


0.5 ±0.2 


1.0 - 


1.5 


0.9 




i^2(1580) 


2- 




1.7 ±0.4 


4.5 - 


6.4 


4.4 




K{U60) 


0- 




forb. 








K*(1410) 




0.8 ±0.4 


4.1 ±0.6 


7.2- 


10.6 


7.3 




iro*(1950) 


0+ 




forb. 








iri(1650) 


1+ 


0.8 ±0.3 


1.7 ±0.6 


not g 


iven 


not given 





Table 4: Comparison of our results for the ratio Rf[%] with previous works. 

8 Summary and Conclusions 

Motivated by the first observation of the radiative decay B 7^^2(1430)7, we have 
investigated rare radiative B decays to orbitally excited i^**-mesons. First, we have pre- 
sented an alternative method of calculating the transition form factors and related decays 
using the QCD sum rules on the light-cone. We restricted our calculations to the leading 
twist-2 operator for the i^'*(892), as in ||22|, and to the asymptotic wave function for the 



other ii'**-mesons. The latter choice is simply based on the fact that using QCD sum 
rules, it is impossible to get rid of the lower-lying states contributions from these higher 
resonances. 

However, nothing is known about the corresponding K**-dec8iy constants, and one has to 
predict them. For that, we have used recent data |^ on semileptonic r — *• (P, S, V, A)^^ de- 



cays to obtain them. For i^'2(1430), we have constrained the corresponding decay constant 
with the recent data H on 5 Kl^lA?,^)-^. We find that if /i^|(i43o) = (140 - 180) MeV , 
a substantial fraction (3.0 — 7.0)% of the inclusive b s'-f branching ratio goes into the 
i('2(1430) channel, in a good agreement with recent CLEO data Our prediction for the 
B — >■ i^'* (892)7 branching fraction yields to (6.0 — 14.0)%, also in good agreement with 
the experimental data As far as decays into higher /T-resonances are concerned, our 

results are in general in much better agreement with |rD| than and [^j, apart from the 



ii'*(1410)-channel where the difference is more significant. For the i^2(1580)-meson, we 
have plotted Br{B — > i^^2(1580)7)/i?r(i? -^s7) as function of the corresponding decay 
constant. 

Finally, it should be noticed that the theoretical uncertainties in our light-cone sum 
rules are the wave functions and the decay constants of the ii'**-mesons. The accuracy of 
our calculation can be substantially improved by taking into account the wave functions of 
twist-3 and twist-4 for the B —>■ J^* (892)7 decay, and going beyond the asymptotic form 
for the other decay modes. To reduce the uncertainties on the i^'**-mesons decay constants, 
one can determine them independently using QCD sum rules for the two-point correlator 
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of the corresponding currents. 
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